Circular solution of two unequal mass particles in post-Minkowski approximation 
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A Fokker action for post-Minkowski approximation with the first post-Newtonian correction is 
introduced in our previous paper, and a solution for the helically symmetric circular orbit is ob- 
tained. We present supplemental results for the circular solution of two unequal mass point-particles. 
Circular solutions for selected mass ratios are found numerically, and analytic formulas in the ex- 
treme mass ratio limit are derived. The leading terms of the analytic formulas agree with the first 
post-Newtonian formulas in this limit. 
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I. INTRODUCTION 

In our previous paper [U (Paper I), a Fokker action 
that describes a time symmetric interaction of point- 
particles is introduced in the framework of a post- 
Minkowski approximation. From two types of Fokker 
action, a parametrization invariant action with a post- 
Newtonian correction and an affinely parametrized ac- 
tion, the equations of motion and expressions for con- 
served energy and angular momentum are derived fol- 
lowing the variational calculation of Ref. Q. We find 
a solution describing a helically symmetric circular or- 
bit in the post-Minkowski approximation (with post- 
Newtonian corrections) that is analogous to the circu- 
lar solution of two charges obtained by Schild for the 
electromagnetic interaction Q. We report here results 
supplementing those of Paper I: numerically computed 
solution sequences for unequal mass particles, and ana- 
lytic formulas in the extreme mass ratio limit. The latter 
results agree with the first post-Newtonian (1PN) formu- 
las; hence a consistency of our model is confirmed in this 
limit. 



II. FORMULAS FOR CIRCULAR SOLUTIONS 

We present a set of formulas governing the helically 
symmetric circular orbits of two point particles, {m, v} 
and {fh,v}, and derive analytic expressions in the ex- 
treme mass ratio limit q := m/fh — > 0. The set of alge- 
braic equations is solved numerically for a fixed binary 
separation to specify each circular orbit. The result for 
the unequal mass binary orbit is presented in Sec. IIIII 
Units of G = c = 1 are used in this report. 



A. Parametrization invariant model with 
post-Newtonian correction 

1. Circular solution 

We first list the result from Paper I for the 
parametrization invariant model with 1PN correction 
terms. The (integrated) equations of motion for parti- 
cles m and fh are written in terms of the velocities, v 



and v, of particles m and m which are related to the 
orbital radius by a :— v/Sl and a := v/Q., 

— m^f 2 v£l = — mfhrf 7Q \F(<p, v, v) 

+ (m + ffi)nFi(tp,v,v,j, ; y)], (1) 

— rfi"/ 2 vfl = —mfh , y"/ 2 ri 2 [F(tp,v,v) 

+(m + fh)QFi(tp,v,v,7,j)]. (2) 

As shown below, {tp, v, v, 7, 7} are not independent. The 
functions F(ip,v,v) — F(tp,v,v) are the post-Minkowski 
terms, while Fi(tp, v, i>,7,7) = Fj((p,v,v,j, 7) is either 
of two alternative 1PN correction terms that agree at 
1PN order: F\ = Fpn(ip,v,v, r y,'y) derived from a non- 
relativistic correction, or Fi — Fgp-^(p, v, v, 7, 7) derived 
from a special relativistically invariant correction. 



F(<p,v,v) 



1 



(ip + vv smtp) 2 



(1 + vv cos ip)v 



x(ipcos<p — v 2 simp) + — y2 )i l P + vvsunp) 

— — [v sin <p(jp + vv sin cp) + (1 + vv cos tp)(v + v cos ip) 



-j—^{¥> + vvsm(p) 2 ]$((p,v,v) \, 
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F PN ((p,v,v,j,j) 



i'SPN (</>,«,", 7; 7) ■= - 



7 2 7 2 (w + v) 3 
1 



1 + -7 v(v + v) 



(4) 



1 



, 3 2 
x< ~T V + 

4 ip + vv smtp 



(77) 5/2 {tp + vv simp) 2 
v sin tp (1 + vv cos ip) (v + v cos tp) 



((p + vv sin tpY 



(5) 



(6) 



The function Q((p, v, v) is defined by 

(1 + vvcosip) 2 - ±(1 - v 2 )(l - v 2 ) 

$(tp,v,v) := ■ — fh . 

ip + vv sin (p 

For the parametrization invariant models, 7 and 7 are de- 
rived from a flat-space normalization of the four- velocity, 

l=(l-v 2 )-K 7 = (l-« 2 )-i (7) 

The retarded angle ip is the positive root of ip 2 — v 2 + 
v 2 + 2vv cos (p. 
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2. Extreme mass ratio limit 

The extreme mass ratio limit q :— m/fh — > is identi- 
cal to the limit v — » with fixed. In the limit v — > 0, we 
may assume that v and fh remain finite. Consequently, 
we have 7 — » 1, cp — > v, and fh — > M , where A/ := m + rn 
is the total mass. With v and O regarded as independent 
variables, Eq. ((TJ is a quadratic equation for OM, whose 
q = form is 



Fi (OM) 2 + F (OM) - v = 0, 
with physical solution 

OM = J-(-F+ ^F 2 + AF lV 



(9) 



The functions F (the post-Minkowski term), F\ — Fp^ 
and F\ — i^sPN (the alternative forms of the 1PN correc- 
tion) for q = become 



F((p,v,v) 



l-3v 2 

v 2 (l — V 2 ) ' 



^spn(^, u,«,7,7) 



(l-v 2 )^f I, 

1 7^ 

4 



(10) 



(11) 



(12) 



where $ has the form $((,2, v, v) = (1 + v 2 )/ (2 w). 

Note that the parametrization invariant post- 
Minkowski model is derived by setting Fj = 0, and there- 
fore OM = v/F. In the q — > limit, this is written 
OM = i; 3 (l-v 2 )/(l-3v 2 ). 



3. Energy and angular momentum formulas 

The conserved energy and angular momentum for the 
parametrization invariant model are written 

E = E PM + ei, and L = L PM + £1, (13) 

where Epm and Lpm are the post-Minkowski terms 

„ 771 fh , 

E PM = - + - (14) 

7 7 



Lpm = 2mm77 v, u), 



(15) 



and ei and £\ are the parametrization invariant 1PN cor- 
rections e\ = epN and £\ — ip^, or those of the special 
relativistically invariant model e\ = espN and l\ = £spn 
given by 

e\ — —£l£i, (f° r b°th e i = e PN and gspn), (16) 



^SPN 



mfh(m + m)0 

77(w + v) 2 
mfh(m + m)0 



(77) 3 / 2 (tp + vvsuiipy 



(17) 
(18) 



In the q — > limit, the conserved energy and angular 
momentum normalized by the mass remain finite. Sub- 
tracting the mass of the heavier particle from the post- 
Minkowski energy, Epm '■= Epm — fh, and taking the 
limit v — > with fh — > M, we have 



(19) 
(20) 
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"« 2 ) 1/2 : 


m 


Lpm 


1 


+ u 2 


mM 


v{l - 


- W 2)l/2' 
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OM, 



OM. 



(21) 
(22) 



^. Solution sequence in q — ► Zimii 

In Paper I, it is proved that the first law of thermody- 
namics that relates the changes in the conserved energy 
and the angular momentum, dE — OeLL, is satisfied by 
binary solutions derived from the parametrization invari- 
ant Fokker action. This relation is used to cross check 
the analytic formula in the q — » limit above as well 
as the numerical solutions shown in the next section by 
calculating dE/dv = QdL/dv, where E := Epm + &i- 

In the parametrization invariant post-Minkowski 
model, the normalized angular velocity of a particle m, 
OM, is defined in an interval < v < l/VS, and it 
becomes infinite at v = l/VS. With the 1PN correc- 
tion Fi — -Fpn, the range of finite OM is approximately 
< v < 0.361598, and with the special relativistic in- 
variant 1PN correction F\ = Fspn, it is < v < 0.36166. 
Newtonian point particles have no innermost stable cir- 
cular orbit (ISCO), but adding a 1PN correction to the 
Newtonian orbit recovers the ISCO that is present in 
the exact theory. In the post-Minkowski framework, 
we find that the existence of an ISCO depends on our 
choice among actions that are equivalent to first post- 
Minkowski order. In particular, we will see that the 
parametrization-invariant action leads to sequences with 
no ISCO even when 1PN terms are included. This is 
plausibly due to the fact that the sequences associated 
with the parametrization-invariant action terminate be- 
fore reaching the angular velocity of an ISCO. In fact, in 
the standard 1PN formalism, an ISCO occurs at an unre- 
alistically high value of angular velocity, namely OM = 
0.544, where the 2PN and 3PN values are OM = 0.124 
and 0.0867 for q — 0, respectively 

Curiously, however, as we note in the next section, se- 
quences associated with the affinely parametrized action 
do have an ISCO, in this case at an unrealistically small 
values of OM. 

Finally, we show that these results of post-Minkowski 
plus 1PN corrections agree with the 1PN formula in the 
q — > limit. In Eq. ©, an expansion of OM in the 
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small v limit becomes f2M = v 3 + 3v 5 + 0(v 7 ) for both 
PN and SPN models, and this is inverted to write v in 
terms of small OM as v = (OA/) 1 / 3 -Q.M +0((CIM) 5 / 3 ) . 
Substituting this into the energy and angular momen- 
tum formulas, the leading two terms agree with the post- 
Newtonian formulas (see e.g. [![) up to the 1PN order 
for the extreme mass ratio q — > 0, 



E 
m 

L 



-(MM) 2 / 3 + -(QM)^ 3 + 0((QM) 2 ), (23) 



mM (flM) 1 / 3 



|(r>M) 2 / 3 + o((nM) i / 3 ) 



(24) 



B. Afflnely parametrized model 

1. Circular solution 

For the affinely parametrized post-Minkowski model, 
analogous forms of Eqs. JT]) and ([2|) are written 

— mj 2 vQ = —mfh'y 2 jft 2 F A (ip,v,v), (25) 

— fhj 2 vCl = —mm r y ; y 2 VL 2 F A ((p,v,v), (26) 
where the function F A (ip, v, v) = F A (ip, v, v) is written 

F A {<p,v,v) := -4 } , ^ 

[ifi + vv sin ip) z 

X /(l + vv cos ip)v((f cos ip — v 2 sin ip) 

+—v(l — v 2 )((p + u-Dsin^) — — [vsinip(tp + vvsinip) 



-(1 + vv cos ip)(v + v cos tp)\ $>((p, v, v) 



(27) 



For the affinely parametrized world line, 7 and 7 satisfy 
- 7 2 (1 - v 2 ) + 4m7 2 7fJ $(<p, u, u) = -1, (28) 

-7 2 (1 - v 2 ) + 4to77 2 u, u) = -1. (29) 

In the limit of q — > (or more directly v — > 0), 
F A (ip,v,v) = (1 - w 2 )/w 2 . From Eq. (J25]) and (28]), we 
have 7 = (1 - w 2 ) 1 / 2 /(l - Av 2 - v 4 ) 1 / 2 , while in Eq. ([2"5]). 
taking v — > and m — * yields 7 — ► 1. As a result we 
have in the extreme mass ratio q — ► 0, 



f2M 



1-w 2 



(30) 



2. Energy and angular momentum formula 

The conserved energy and angular momentum for the 
affinely parametrized model are written 

vn TYi 

E = h — + 4<mfh'yjQ,$((p,v,v), (31) 

7 7 




FIG. 1: Angular velocity, in dimensionless form QM, is plot- 
ted against the velocity of the lighter particle for 3 mass ratios 
and the q — » limit for the parametrization invariant model 
with SPN correction. Curves of the analytic solution for q — > 
and that of q — 0.001 overlap each other in the plot. 
1 
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FIG. 2: Energy, in dimensionless form E/m, where E — E 
rh, is plotted against QM for the same models as in Fig[T] 



where the form of $(<p, v,v) is the same as that of the 
parametrization invariant model ([6]). Using Eq. (|28p and 
([29|) . the energy can be rewritten 

lml , o . 1 fh 1 . 0s 
E = -- + -m 7 1 - v 2 ) + r-r + ^7(1 - « 2 )- (33) 
z 7 z z 7 z 

In the g — -> limit, the energy without the rest mass 
of the heavier particle, E :— E — fh, and the angular 
momentum become 



E 
m 



(1 - 3v 2 ) 



{1 - v 2 )(l - 4v 2 - v*)]V 2 



L 1 + v 2 / 1-v 



mM v V 1 — 4w 2 — v 4 



,,2 \ V2 



(34) 
(35) 



L = 2mm77 $(95, u, u), 



(32) 
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FIG. 3: Same as Fig. [TJ but for the affinely parametrized 
model. In the q = 1 case the ISCO occurs at v ~0.184. 
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gular momentum have a simultaneous minima at v = 
v /(l + 2 4 /3_25/3)/3 pa 0.339136, which corresponds to 
flM « 0.0440743. 

III. NUMERICAL SOLUTIONS FOR THE 
UNEQUAL MASS CIRCULAR ORBIT 

A circular solution is calculated from algebraic equa- 
tions given in Sec lII A Tl for the parametrization invariant 
model, and III B Tl for the affinely parametrized model. It 
turned out that a convenient way to find a solution is (1) 
fix the ratio of velocities v/v and determine the corre- 
sponding mass ratio from the equations of motion, then 
(2) change the velocity ratio to adjust the value of the 
mass ratio to a fixed value (using the bisection method, 
for example). 

In Figs. [JJ and the plots of the parametrization in- 
variant model with SPN correction terms are presented 
for three mass ratios, q = 1.0, 0.1, and 0.001. Plots for 
the case with PN correction terms are not shown here, 
and they are qualitatively the same as SPN cases. The 
analytic solution in the q — > limit is also plotted and it 
overlaps with the q = 0.001 line in the plots. 

In Figs. [3] and |4j the plots of the affinely parametrized 
model are presented for the same mass ratios as above. 
The solutions of the affinely parametrized model are 
markedly different; for any mass ratio q € [0, 1], we found 
a simultaneous minima in the energy and angular mo- 
mentum which corresponds to the ISCO. 
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QM 

FIG. 4: Angular momentum, in dimensionless form J/(mM), 
is plotted against QM for the same models as in Fig. [3] Min- 
ima of each curve corresponds to the ISCO. 

3. Solution sequence in q — > limit 

The first law SE = flSL is also satisfied for the affinely 
parametrized model, and hence one can cross check for- 
mulas in the q — > limit using the relation dE/dv = 
CldL/dv. Although the normalized angular velocity of a 
particle m, VIM, is finite in an interval v £ [0,1), the 
redshift factor 7 as well as conserved quantities E and L 
become infinite at v = \J \7I — 2 « 0.485868, which cor- 
responds to QM = (VE- 2) 3/2 /(3 - V5) w 0.150142. 

In this interval, v £ [0, \J \/5 — 2), the energy and an- 



IV. DISCUSSION 

Agreement between the energy and angular momen- 
tum formulas of the 1PN circular solution, and those 
of the parametrization invariant post-Minkowski model 
with post- Newtonian correction, is exhibited only for the 
extreme mass ratio limit in this report. For an arbi- 
trary mass ratio one needs to expand the retarded an- 
gle tp to the next order in the velocities, v and v, as 
if w (v + v)(l — vv/2), and the rest of the calculation 
closely parallels that of the q = case. 
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